Two dimensional electron gas driven by strong alternating electric field 
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2D Fermi liquid driven by uniform ac-electric field at zero temperature may remain in quantum 
coherent non-equilibrium state. We develop a quasistatic approximation for strong and slow ac- 
fields and solve the problem of driven disordered 2D electron gas in high non-overlapping Landau 
levels. The broadening of Landau level has the Lorentz form and is inversely proportional to the 
amplitude of ac-field. In the absence of electron-phonon interaction the electron distribution function 
is constant within the last Landau level and the diagonal dc-conductivity is zero. For weak electron- 
phonon interaction the rfc-conductivity is anisotropic. A kinetic transition from the phonon cooling 
to the phonon heating is predicted. 

PACS numbers: 72.10.-d, 73.40.-c, 78.67.-n, 73.43.-f 
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I. INTRODUCTION 

Experience with microwave ovens teaches us that elec- 
tron systems subjected to an alternating (ac) electric 
field are heating up. This paper explores examples to 
the contrary. A closed system remains in a quantum co- 
herent state even if driven by ac-field. This driven state 
is equivalent to the static ground state provided the sys- 
tem possesses the translational invariance, the parabolic 
band dispersion and instantaneous electron-electron in- 
teraction. Boundaries and disorder preserve the coher- 
ence but profoundly transform the driven state. It is 
the electron-phonon interaction that destroys the coher- 
ence and eventually heats up the electron system. At low 
temperature the electron-phonon interaction is especially 
weak and the quantum coherence of a driven system can 
be observed. In three dimensions the ac-electric field is 
confined near the surface but in two dimensions (2D) it 
penetrates and drives electron system into a new state. 

Recently in high mobility GaAs heterostructures un- 
der microwave irradiation in perpendicular magnetic field 
new oscillations of conductivity and new zero-resistance 
states has been discoverccl^-^. These oscillations depend 
on the ratio of the ac-frequency lj to the cyclotron fre- 
quency ujh- Ryzhii has predicted a phenomenon of the 
absolute negative conductivity in two dimensional elec- 
tron gas (2DEG) driven by the microwave field^ and, re- 
cently, his theory has been confirmed in several models^ 
of disordered electron system. 

Proper description of the absolute negative conduc- 
tivity requires to account for a current necessitated by 
electron relaxation from the exited Landau levels pre- 
dominantly by emitting phonons^. This compensating 
current amounts to a positive conductivity. In GaAs the 
phonon effect is weaker than the Ryzhii effect of disorder 
and the theory neglects phonons in the limit ojh <C Eph^ 
where Eph is the characteristic phonon energy (5 Kelvins 
in GaAs). 

The problem of disordered 2DEG in magnetic fields 
that quantize the electron motion is a special one because 
it is the disorder that lifts the Landau level degeneracy 
and imparts a velocity to electrons, and it is the disorder 



that is coupled to the ac-field. The quantum limit cor- 
responds to a weak disorder: loht 2> 1, where r is the 
elastic scattering time. In high Landau levels the exact 
self-consistent Born approximation theory predicts the 
semi-circle density of states and longitudinal resistivity 
proportional to the magnetic field''. 

In this paper I find modifications in this picture due 
to microwaves in the limit of slow and strong ac-field E: 
Lo <^ luh and uj^ ^ {eEvpYT, where vp is the Fermi ve- 
locity, called a quasistatic approximation. In strong fields 
an electron drifts with large velocity 'over the disorder' 
and experience predominantly forward scattering^. Ac- 
cordingly the longitudinal rfc-conductivity of the driven 
2DEG is strongly suppressed and is zero in the absence of 
the electron-phonon interaction. The semi-circle density 
of states of the static 2DEG transforms into the Lorentz 
spectral density of the driven 2DEG with the broadening 
of Landau levels being inversely proportional to E. 

In the first two sections we overview a Fermi liquid 
driven by ac-field. Here we set a stage by examining 
examples from the usual Fermi liquid. In particular a 
quasistatic approximation and a special driven state of 
ID Fermi gas with boundaries that possesses a smooth 
electron distribution function are introduced. In Sec. Ill 
we find a kinetic transition from the phonon cooling to 
the phonon heating. In the last two sections we solve 
the problem of the ac-drivcn disordered 2DEG in high 
non-overlapping Landau levels with Sec. V containing 
the main results. 



II. FERMI LIQUID IN ^C-ELECTRIC FIELD 

In this section we overview the theory of Fermi liquid 
driven by uniform ac-electric field. We determine the 
quasiparticle dispersion, the linear response function and 
effects due to interactions with phonons and impurities. 

Uniform in-plane electric field E{t) can be described 
in the gauge with zero electrostatic and uniform vector 
potential A{t) = —c J E{t)dt assumed to be limited and, 
therefore, with no arbitrary gauge degrees of freedom. 
We write the Hamiltonian of 2D Fermi liquid H{t) = 
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Ho[i'(t)] as: 

and put ?i = 1 throughout this paper. Classically a free 
electron oscillates in ac-field with the velocity: V{t) — 
—eA{t)/mc. The wave-function ^'(t) of a closed system 
is a solution of the Schroedinger equation: 

z^^^it)^ limit). (2) 

If the Hamiltonian does not depend on time then its 
eigenvalue - the energy - is conserved. The momentum 
operator commutes with the one-particle Hamiltonian 
Hq [A] and the momentum of electron p is a good quan- 
tum number. The solution of one-particle Schroedinger 
equation with momentum p yields: 

^,(t,.-o.-^exp(^-y_(p+-A(t)) dty (3) 

where S is the area of sample. If A{t) is a periodic 
function of time with the period 27r/cj then apart from 
the multiplicative factor exp(— ief) the electron states Q 
- known as Floquet eigenstates - are periodic in time 
whereas e — 121X1 mod (0 < e < is a quasienergy. 
We expand the Hamiltonian iQJ in powers of A{t): 

H{t) =Ho + H"'' - -Ait)^ + -^A^t)No, (4) 
c Zmc 

where Hq — Hq [0] and Hq + ™* is the static Hamilto- 
nian, Jo — {ie/mc)J2i'^i and iVo = X^i 1 ^''^ ^^'^ 
tal electron current and number operators. All these 
three operators commutes with each other. Therefore 
the many-body wave function can be chosen to be the 
eigenfunction of all these three operators with the total 
energy E, the total number of electrons and the total 
momentum P and can be expanded in the basis of the 
momentum eigenfunctions: ipp{r) = exp(ipr)/\/!S,: 

N 

*o = e-'^*^ a(pi..p^)J]Vp-(rl), (5) 

P1--PN *=1 

where a(pi..pV) are complex- valued coefficients that arc 
non-zero if ^^Pi — P- The wave function: 

N 

VI/ (i) = e-'^* n (6) 

P1--PN i=l 

where il)p^{t,ri) is given by Eq.Q, is the solution of 
the time dependent Schroedinger equation |(5J with the 
Hamiltonian It is the eigenfunction of the to- 

tal electron current and number operators: Jo^(i) = 
eP/mc and iVo*(t) = A^*(i). The solution ® 



requires the same conditions as the Kohn theorem^: 
the electron band dispersion is strictly quadratic, the 
electron-electron interaction is instantaneous and the sys- 
tem is translationally invariant without boundaries or 
bulk disorder. Weak interactions violating these condi- 
tions - boundary, disorder or the electron-phonon inter- 
action - can be treated in perturbation thcory'\ 

Canonical transformation: U'^ {t)[id / dt — H{t)\(j{t) ~ 
id/dt— Hq — if'"*, of the second quantized electron op- 
erators: il'ip) — U{t,p)x{p), where 

(7) 

eliminates the ac-electric field from the Schroedinger 
equation. Essentially it is a transformation from the lab- 
oratory frame to the center of mass frame oscillating with 
electrons^. If a stimulus is canonically transformed then 
the response is the same as in the static system. 

The electron density operator p{t,r) — ip^r^ipir) in 
the driven 2DEG can be expressed in terms of the static 
density operator by using the wave function 0: 

P{t, q) = pM cxp [iq j V{t)dt^ (8) 

where /3o{g) = X]p + O) is the static density 

operator. Therefore, a generic instantaneous transla- 
tionally invariant electron-electron interaction: iJ"** — 
^'YllqPo{q)U{<t}PG{~'q)^ is invariant under the canonical 

transformation 0: i/"* = \Y.qP{i^q)U{q)p{t,~q)- 

In the static Fermi liquid the quasiparticle dispersion 
reads: ^(p) — w_f(|p| — pF)r^ where pp is the Fermi mo- 
mentum related to electron density. In the driven Fermi 
liquid the Floquet theorem and the momentum conser- 
vation allows for two quasiparticles from the inside of the 
Fermi circle to scatter into empty states outside the Fermi 
circle with the total energy increment being an integer- 
multiple of Lj. But the existence of canonical transforma- 
tion Q shows that the matrix elements of such events 
are zero and that the quasienergy of quasiparticle is the 
same: ^(p) = u_f(|p| —pp). In Landau theory the relax- 
ation of quasiparticles with momenta close to pi? is due to 
the electron-electron interaction and in both static and 
driven Fermi liquids Tin ^ (p ~ Pf)~^ ■ 

Incoherent driven state. The static ground state is 
transformed into a driven state during adiabatic switch- 
ing of the ac-field at zero temperature. This driven state 
is no longer the unique lowest energy state because the 
energy is no longer conserved. In the absence of the 
electron-phonon interaction the driven state is a quan- 
tum coherent state but weak electron-phonon interaction 
directs the driven system into an incoherent state char- 
acterized by an effective electron temperature T*. The 
electron-electron scattering probabilities are the same in 
both driven and static Fermi liquids. Therefore the solu- 
tion of the Boltzmann kinetic equation with the electron- 
electron collision integral yields the Fermi-Dirac distribu- 
tion function /(e/T*) that depends on the eigenvalue of 
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the static Hamiltonian e and some T* conserved by the 
electron-electron interaction, despite the bulk being kept 
at r = 0. T* is determined from the balance of two 
energy flows: the 'heating' disorder scattering and the 
'cooling' phonon emission. 

Coherent quantum state of driven system with 
boundary emerges in the absence of electron-phonon 
interaction. Consider ID non-interacting electron gas 
driven by parallel ac-field and confined to a box of length 
L with two boundaries of zero transparency. Electron 
states near the Fermi points ±pf ^ 1/L are the sum 
of wave functions Q with left and right moving compo- 
nents: 



E 

p>0 



,ip(x~X(t)) 



B{p)e 



-ip(x~X(t)) 



(9) 

where X{t) = J V{t)dt and the ac-field is periodic: 
V{t) = Vsm{ujt). Floquet theorem states that momenta 
in Eq.ljS)) are quantized: = Pq + 2mnuj, where pq is 
the middle momentum and n is integer. Correspond- 
ing amplitudes A„ — A{pn) and i?„ — B(pn) are re- 
lated as An = (— 1)"+-'^_B„. Two boundary conditions: 
^U=-L/2(i) = and ^\x=L/2{i) = 0, can be rewritten 
as: MnkAk ~ 0, with the matrix 



Mnk = z"-V„_fc (^V) Sin + , (10) 

where Jn{x) is the Bessel function, po is found from 
the self-consistency equation detM(po) = 0. We solve 
numerically for An and find a broad electron distribu- 
tion over the energy range: Se ~ nmasW ~ Vpp- A 
Slater determinant constructed from all wave functions 
with \pq\ < pf is the quantum coherent state that turns 
into the static Fermi gas at ^ = and that describes a 
smooth electron distribution in the momentum space. 

Linear response function of charged Fermi gas is 
the response to the scalar potential </)(t, f). If the ac-field 
velocity is periodic: V{t) = V sm.{ujt), then the response 
function: 



PkL 



{p{t,q)) ^X(t,t',q)^{t',q), 



(11) 



is also periodic: (/5(r2 + ncj, q)) = x„(ri, q)4>(Q., q), where 
il is the frequency of field (j>{t,f). We use the Keldysh 
method for non-equilibrium systemsii to find the gener- 
alized Kubo formula: 



27r Jo 



intuT 



dT 



{[p{q,T~t/2) p{~q,T + t/2)l 



dt 



jnt 



(12) 



where the brackets stand for the matrix element over 
the driven state with electrons filling the Fermi circle of 
radius Pi? and the density operators p are given by Eq. ^ . 
After taking the integral with respect to T and t we find: 



A;— — oo 



1 



1 



P{n + qy2m)-P{n-qy2m) \ 
I- I —/ ' (13) 



where P{x) = \/ v\(p — (a; -)- (fc + n/2)LoY. The static 
response function for lo <^ vpq is found as: 



X2n{q) = —Jl[Vq/^ 



(14) 



Zeroes of this response is an analog of the phenomenon 
of the coherent destruction of timneling^^. 



III. QUASISTATIC APPROXIMATION 

In this section we introduce a quasistatic approxima- 
tion and we discuss the role of disorder and electron- 
phonon interaction. 

For slow electric fields it is convenient to take the 
Fourier transform of Ea. l|ll|) with respect to the 'fast' 
time t — t': {p{fl,q,t)) ~ xi^, O)- For strong 

electric fields: to <^ PfV, a large number of harmonics 
contributes to the sum 1)1311 . In this limit the response 
function 1)13(1 has an asymptote: 



Ztt \ q^/m 



(15) 



where P{x) = ^ vj,q^ ~ {x + V{t)qj . Eqs.^MB be- 
come identical in the limit uj pfV <C cf and q -^pf- 
The Green functions in the laboratory and in the 
center of mass frames are defined as: Giab(t,t' ,p} = 

-i{Tti;{p,t)ip''{p,t')) and 



G{t,t',p)^-t{Tx{p,t)xHp,t')), 



(16) 



where T is the chronological ordering operator. These 
two Green functions are related by the canonical trans- 
formation ^ . The Green function in the center of mass 
frame ()16(l depends on t — t' and is determined by the 
Hamiltonian Hq + iJ*"* of the static Fermi liquid: 



Gie,p) 



e-^(p) +iv{p- Pf)\p- Pf\' 



(17) 



where a, rj are the amplitudes of the pole renormalization 
and relaxation^°. In the Fermi gas a = 1 and r] — > -1-0. 

Disordered Fermi gas (non-interacting) can be de- 
scribed by the model of weak short ranged impurity po- 
tential treated in the Born approximation. The disorder 
Hamiltonian reads: H™v — J U{'r)p{t,r) d^f, wheie the 
random Gaussian disorder potential has the correlation 
function: ([/(0)[/(r)) = 5{r)/2-KVFT. r is the mean elas- 
tic scattering time and z/;- = m/27r is the 2D density of 
states. Disorder breaks down the translational symmetry 
and matrix elements of the disorder Hamiltonian become 
time dependent in the center of mass frame. On the other 
hand the Green function averaged over the disorder is 
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translationally invarianlr^ both in the static and driven 
Fermi gases. In the standard diagrammatic approach^" 
that neglects diagrams with crossing impurity hnes (small 
as l/epT) and in the limit uj <^ \t — <C cf we find 
the electron self-energy: 



S(t-t',p) 



JV{t)(p'-p)(t-t') G{t ~ t',p') (fp' 
2-KVFT (27r)2 ' 



(18) 

depending on the time i as a parameter. Taking the 
Fourier transform of Eq.lQHJ with respect to t — t' and 
using the definition of the self-energy: G"^ (e) = e— ^(p) — 
E(e,p), where ^{p) = p^ /2m — /i, we find: 
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e+-e'-S(e+,p') 



de', (19) 



where e+ ^ e + V {t){p' — p) . The solution of the Dyson 
equation ifT?^ yields: T,{€,t,p) — —isgn{€~V{t)p)/2T 
(neglecting the time dependent real part) and the Green 
function in the center of mass frame: 



Gie,p) 



a?) + I sgn( e-V{t)p] /2t 



(20) 



The corresponding Green function in the laboratory 
frame is: G,afc(e,p) = (e + V{t)p- ^{p) + i/2T sgn(e))-i 

One Green function 12()|l describes two distinct situa- 
tions. In the absence of electron-phonon interaction the 
electron distribution function f{p) = 1 if e — V{t)p < 
and zero otherwise. In the center of mass frame this 
Fermi circle oscillates slowly. On the other hand when 
electrons could quickly lose the energy into the bulk 
the driven state is incoherent with f{p)~l inside the 
stationary Fermi circle. Here scattering between few 
electron states close to the Fermi circle with the re- 
versed imaginary parts: e — V{t)p < and p > pi? or 
e — V{t)p> and p < pp, creates the Joule heating. 

Phonon cooling is conveniently described by the 
electron-phonon Hamiltonian in the deformation poten- 
tial approximation (Froelich model): 



(21) 



where O is the isotropic deformation potential and p{f, t) 
is given by Eq.©. Electron self-energy is found in the 
second order of the electron-phonon interaction H21|l : 

+ (22) 

where p is the crystal density, uj = lu + V{t)q, q is the 
3D phonon wave vector. There is in- plane momentum 
conservation: g|| =p'—p. The velocity of free electrons in 
ac-field is V{t). The phonon propagation for zero lattice 
temperature is retarded: 



D{LU,q) 



w2 - (sg)2 + iO 



(23) 



We transform the measure of integration in Ea. H22|l . Let 
(j){q) be the angle between the momentum p' and the mo- 
mentum p. It can be found from the equation: p'^ +p'^ — 
2pp' cos (t){q) = g^. We substitute (f){q) in the measure of 

integration (Pp' = p'dp'd(l)(q) and find the transformed 
measure: q\\p' dp' dq\\/ ^ p']pq^^ — 911/4— (p^ —p''^Y/A. In 
the hmit p^ — p'^ <C PFq\\ Ea. (P3 can be rewritten as 



^ vps" j G{e + ijj,i')D{Cj,q) dujdq\\dqzd(,' ^^^^ 



El, 



PFJl-qfJ'ipj, 



(2^)2 



where Eph = \/ ps^ /Q^ is the characteristic phonon en- 
ergy. Imaginary part of the self-energy gives the prob- 
ability of the ac-field assisted phonon emission. Multi- 
plying it by the energy transfer t^(t)qj| — sq we find the 
cooling rate. Electron-electron interaction establishes the 
Fermi-Dirac electron distribution function: f{e/T*) = 
l/(exp(e/T*) -I- 1) with the electron temperature T* de- 
spite the bulk being kept at zero temperature. 

Following the standard method^*' we find the following 
results. As long as V{t) < s there is no energy absorption 
from the ac-field into the electron system associated with 
the phonon emission. Using the identity: 



/(f 



1-/ 



T 



de = LuN 



where N{uj/T) = l/{exp{bj/T) - 1) is the Bose-Einstein 
distribution function, we find cooling rate per unit area: 



dt 



2pFEf 



ph 



iSsq?N(%l]'^^, (26) 



{2n) 



where 6s = s — V (t) cos{9) cos{(f>) , q^ = qsinO, <7|| — 

qcosO and cj) is the angle between vectors V{t) and 
Evaluating the integral we find the cooling rate: 



dE 
'dt 



12C(5)sz.|, r5 / K(fc) , (3 - fc2)E(fc) 



2^'pfE'p, 



l-F 



(1 



(27) 

where k — V{t)/s. The cooling rate diverges as V{t) s. 

At V{t) > s a channel for heating is opened. Electrons 
in a narrow area in momentum space emit phonons and 
absorb the ac-field energy. This part of the heating rate: 



dE _ 64p|, 
It ~ 90ttE^ 



{V{t)-sfp%eiV{t)-s), (28) 



ph 

where 0{x) is the step-function, is additive to the coohng 
rate (|27|l from the rest of electron states. Heating H28(l 
does not depend on T* as long as it is small. And the bal- 
ance of energy flows Eas. (|27l28|l makes the electron tem- 
perature to grow rapidly: T* ~ \V {t) - sl^/'^Vp^^ s~^/^pf- 
Thus there is a kinetic transition from the fast cooling 
Ea. (|27|) to the growing heating Ea. (|28|l at the critical 
value of ac-field: V{t) = s. 
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IV. ELECTRON IN CROSSED MAGNETIC 
AND ^C-ELECTRIC FIELDS 

In this section we find the wave function, density and 
current operators of free electrons in crossed magnetic 
and periodic ac-electric fields. 

Consider an electron confined to a plane in magnetic 
field H perpendicular to the plane and driven by in-plane 
electric field E{t) — E cos{Lut) that is linearly polarized 
in X direction. In Landau gauge: = —Hy , and in the 
vector potential gauge for the uniform ac-field: A^lt) = 
—cEsm(Lut)/uj, the Hamiltonian reads: 



1 

2m 



-^^ + li-Hy + A,it))) (29) 



dx 



2m dy"^ 



In high Landau level N there are two distinct lengths: the 
cyclotr on radiu s Rc = V2NIh and the magnetic length 
Ih ~ ^Jhc/eH. We use the magnetic units where e — c, 
H = 1 and Ih = 1. A characteristic energy of the ac-field 
is defined as: £ = eERc- The electron wave function in 
Landau level N and in the state of the gauge index p:^^ 



ipNp{t,r) 
exp{iX{t){y - 



CKp{ipx) 



p-Y{t))-{y-p 



HN{y-p-Y{t)) 



r(^))V2 + ^$)(30) 



is the solution of the time dependent Schroedinger equa- 
tion with the Hamiltonian H29|l . X{t) and Y{t) are the 
X, y-coordinates of the electron classical trajectory - an 
ellipsoid - in the crossed electric and magnetic fields: 



Xit) 



Eloh cos{ujt) 



Y{t) 



where ujh = eH/mc, loj^ = cj^ ± Ltrfj and: 



sin(2wi) - ujut \ N + - 



(31) 



(32) 



is the wave function phase. The retarded bare elec- 
tron Green function is diagonal in the indices N and p: 
gR{e) = l/(e - Nujh + iS), where < S 0. 

Using the wave function (|30() we find the density oper- 
ator p{t, r) = '(/'^(t, r)ip{t, r) of the driven 2DEG as: 

P{t^ = Mq) exp {-iqyX{t) + iqxY(t)) , (33) 

where the density operator of the stationary 2DEG reads: 

nn' p 

(34) 

The vertex Vnn' [q) for electron scattering from the Lan- 
dau level n to the Landau level n' with transferred mo- 
mentum q is defined as: 



Vnn'{q) = ( -ZT 



l\\ ±1/2 / I ■ \ \n—n'\ / 2 

''N ' qy±iqx\ ' |n-n'| ((£_ 

min(nn') I 2 



V2 



where L^^ix) is the Laguerre polynomial. The two signs 
± in Ea. (|35|l corresponds to the two cases n > n' and 
n < n' . 

Instantaneous electron-electron interaction is invariant 
under the transformation H33|l of the density operators 
from the stationary to the driven form for any time de- 
pendence of the ac-field. For instance the Laughlin wave 
function is as good for driven 2DEG as for the stationary 
one. Variation of the Hamiltonian H29() with respect to 
the vector potential 5A[r) gives the current operator: 



= H + {iujH cos(wt) -I- oj sin(wt)) £/2V7Vwi (36) 

where the operator H lowers the Landau level state 
H|7V,p) = ^/N\N -l,p). P{t,r) is a Hermitean con- 
jugated operator to J{t,f). In Sec. V we neglect the 
second term in the r.h.s. of Ea. (|36|l as being small as 
in the high Landau levels. We will study effects of 
ac-field in the leading order of l/N'^. 



V. DISORDERED 2DEG IN THE QUANTUM 
LIMIT 

In perpendicular magnetic field that quantizes the elec- 
tron motion, electron states are degenerate and it is the 
disorder that lifts this degeneracy and imparts a veloc- 
ity to electrons. Because disorder couples to ac-field the 
problem of Landau level broadening and dc-conductivity 
has to be reexamined. In this section we solve this prob- 
lem in high non-overlapping Landau levels and for short- 
range correlated Gaussian disorder. 

The Boltzmann kinetic equation is useless for calculat- 
ing the conductivity in the quantum limit. The current 
is the sum: (j) — e^^ViSfi, over electron states i and 
a variation of the electron distribution function in weak 
dc-electric field is 5fi ViE. Because ui is random in the 
quantum limit it follows that (v) = and {6f) = 0. In- 
stead we adopt a diagrammatic approach that calculates 
an average of the diffusion coefhcient Di ^ vf. 

We assume that the transport of 2DEG can be prop- 
erly described in terms of the non-interacting electrons. 
In the framework of Landau theory Eliashberg has proved 
that the conductivity of Fermi liquid is the conduc- 
tivity of non-interacting quasiparticles near the Fermi 
surfacaiSi. In magnetic fields the de Haas-van Alphen 
oscillations of the magnetic moment has also been de- 
scribed by the non-interacting quasiparticlea^^. 

The impurity scattering dominates the transport at 
low temperature. The Hamiltonian of electron moving in 
a random impurity potential U{r) reads: 



(35) 



(37) 



where we omit the spin index of the creation and annihi- 
lation operators: V'^ and ^ due to the vanishing exchange 
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interaction in high Landau levels. The probability distri- 
bution functional of the impurity potential U (f) is as- 
sumed to be Gaussian with the correlation function: 



{U{r)U{r')) = 2'Kul]j5 r'^ , where u = loh/t, 



where r is the mean elastic scattering time. 
We expand the electron Green function: 

G{t,t'-r,r') = ~i{Ti;(t,r)i;\t',?)), 



(38) 



(39) 



into a series over the impurity potential U (r) and then 
average it over the disorder, using the standard diagram- 
matic method of crossed diagrams^". In Landau level 
the coordinate dependence of the average Green function 
G{t,t';f,r') — GN{t,t')^]y{f,r') is given by the univer- 
sal phase factor $Ar(r, r') which is conveniently trans- 
fered onto the vertex elements Ea. (|34l35|l . GN{t,t') is 
the sum of diagrams consisting of bare electron 5Ar(e) 
and impurity 138|) lines incident in vertices (see Ea. (|35|l ): 
Vnn'iq)exp{iqx{p + qy/2) - q"^ / 4)6p' ^p+g^ , that describe 
a process of electron scattering off the impurity poten- 
tial U{q) from the state {n,p) into a state {n',p') with 
transfered momentum q. 

In the absence of ac-field the density of states in a 
high Landau level follows the semi-circle law: p{e) — 
\/4u — e^/2TTl'jjU, where e is zero in the center of the 
Landau level. At T = the electron distribution 
function /(e) is given by the step function with the 
Fermi energy being determined by the electron density 
n = {N + v)/2-kI'\j, where v is the filling factor of the 
last A^th Landau level. The longitudinal conductivity: 
r^xx = ne'^Nup'^{€F)r^ is comparable to the Hall con- 
ductivity a^y = ecn/H. 

Consider 2DEG driven by strong and slow linearly- 
polarized ac-elcctric field: £ — CxS cos{ujt), where <C 
£ <^ loh and u <^ u/£. We neglect Landau level mixing 
small as ^/u/loh and impurity lines crossings small as 
and omit the index N from the notations. The self- 
energy is given by the sum of all non-crossing diagrams: 

Y.{t,t')^u{t,t')G{t,t'), S(e,i)= [^u{n,t)G{e + n,t) 

J 2tt 

(40) 

The impurity line includes the correlation function 138|) 
and the two adjacent vertices Ea. (|33l34l35|l : 



{2nY 



(41) 



where 5r{tt') — R{t) — R{t') is the differential along the 
classical electron trajectory (|31ll . The integration with 
respect to the transferred momentum q gives: 

u{t,t') = u [L% ((5r2(tt')/2)]'exp(-<5r2(tt')/2) , (42) 

where U^(x) is the Laguerre polynomc. In the limit: 
TV > 1 and u)\t - t'\ <C 1, Eq.(|^ has an asymptote: 



where c{t) = y cos^ ujt + {lu/ujh)'^ sin"^ ujt. The Fourier 
transform of the impurity line with respect to the time 
t — t' reads: 



'1- 



4f2c2(t) 



(44) 



where A{t) = u/-K£c{t) is the characteristic broadening 
of the Landau level, K(fc) is the Jackobi elliptic function 
of the first kind. u{Vt,t) = for \n\ > 2£c{t). The non- 
crossing approximation is valid provided 5r{tt') <^ Ih or 
£c[t)\t - t'\ < %/lV. In the limit 1 < £c{t)\t - t'\ we find 
approximately: u{t,t') = A{t)/\t — t'\, where oscillations 
with the frequency e„p = 2£c{t) are neglected as resulting 
in small corrections of the order of -Jul£c{t). 

For a free electron in random potential the Green 
function G{t,t') can be chosen retarded with the prop- 
erty G^{t,t') ^ a.t t < t'. Therefore we substitute 
\t -t'\ ^ t - t' in u{t,t') = A{t)/\t - t'\ and take the 
Fourier transform of the self-energy H40|l with respect to 
the 'fast' time t — t': 



de 



iAit)G%), 



(45) 



where the 'slow' time t serves as a parameter and we omit 
it from the notations. Solution of the equation G~^(e) — 
e — S](e) together with the Dyson equation H45|) yields 
the Green function having an asymptote: G^{e) = 1/e 
as |e| — » oo, implicitly as: 



1 



G«(e 



iA{t) log 



£c{t) 



A{t)-zG],'ie) 



(46) 



For small energies e ~ A{t) <^ £c{t) we find G^^(e) ~ e+ 
iA{t) log {£c{t)/A{t)). Advanced Green function is com- 
plex conjugated to the retarded one: G^(e) = {G^{e))* . 

In driven 2DEG p(e) = ImG^(e)/ TT is the spectral 
density of electron states - an analog of the density of 
states in the stationary 2DEG. The spectral density 146|) 
has an asymptote: ImG^(e) = A{t)\og{£c{t)/e) /e^, as 
A(i) ^ |e| <^ £c{t). At small e the spectral density has 
approximately the Lorentz form: ImG^(e) = AL/{e^ + 
A^L^), where L = log (f c(<)/A(i)). At high frequen- 
cies |e| > 2£c{t) we find a crossover tail ImG^(e) ^ 
A^/e^dej — 2£c{t)) towards the self-consistent Born ap- 
proximation solution p(e) = for |e| > 2^/u. Using 
Ea. (|45|) we prove that the spectral density is normalized: 



J ImG^ (e) de 



(47) 



u{t,t') =uJ^{£c{t){t-t'j), 



(43) 



We substitute G^de = dG^/{Gj^ + iA) into the inte- 
grand 2i\mG^de = G^de—G^de and deform the contour 
of integration to the real axis: dG~^ /{G^ + iA) — 
dG^^HG-/ - iA) = -2i J dGA/{G^ + A^) = -2iTr. 
Retarded or advanced Green function alone could not 
describe many electron systems and the electron distri- 
bution function F{e) has to be determined as well. We 
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l+l 



a 



n 



n 



n+1/ 




in Fig.l (where uj is no longer the ac- frequency): 



FIG. 1: Two diagrams: a and b, that give the impurity con- 
ductivity in strong magnetic fields 



derive F{e) from the average Keldysh Green functionii. 
We sum up non-crossing diagrams using the Keldysh 
methodii and find the Dyson equation for F{e): 

[F{e) - f{e)] ImG«(6) = G^{e)G^{e) J ^ 

u{n) [F{e + /(e)] ImG^(e + n), (48) 

where Green functions are given by Ea. H4t)|) . At T = 
the initial electron distribution function /(e) is the Fermi- 
Dirac step function. We prove the identity 

—u{n)ImG^{e + n)= ImG^ie), (49) 
27r 

by evaluating the imaginary parts of the two equations: 
i) G(e) = l/(e- I](e)) and ii) the Dyson Eq.^. From 
i) we find that ImG(e) = G^(e)G^(e)ImE(e). Then from 
the Dyson equation (|40|l the identity (|49|1 follows. A con- 
sequence of this identity is that the solution of the Dyson 
equation H48|) does not depend on the initial distribution 
function /(e). Solving Ea. l|48() for F{e) and using the 
normalization of the total number of electrons we find 
that F{e) = v, where ly is the filling factor of the iVth 
Landau level. This solution interpolates between F = 1 
and F = for TV — 1 and -|- 1 Landau levels making 
connection in the gap regions where p(e) = 0. The con- 
stant distribution function is an analog of the classical 
energy diffusion for the non-stationary scattering.^'' 

We use the Keldysh method^ ^ to find the diagonal 
conductivity. It is explicitly a real function: axxiS^) — 
axxi—^)*, with the reactive and the inductive parts: 
Z{rt) — R{il'^) +inL{il'^) , and we disregard the inductive 
part here. For the reactive part we write: 



,o^ 2 2 r de Fie)-F{e-n) ^ 



where the conductivity of the states with energy e is: 

oc 
nn' — 1 

ImG„_i^„'_i(e) + ImG„„'(e)ImG„_i,„'_i(e - f^))(51) 

that includes inter Landau level transitions. The disorder 
average (...) is expanded into a series of diagrams. In the 
leading order the conductivity is the sum of two diagrams 



a^xi^) ± (^yy{^) = '2e^N / / ImG(e 



ImG(e) 



P{t,uj) \ F{e + uj-n)- F(e) deduj 



Sit,Lu) 



{2n) 



(52) 



Functions P(t, uj) and S{t, to) are given by the sum of the 
two impurity lines on diagrams a) and b) in Fig.l. P(t, uj) 
corresponds to diagrams with the two different J and 
vertices (shown in Fig.l) whereas S{t,Lo) corresponds to 
the diagrams with J ~ J or — pair of vertices (not 
shown): P{t,t-t') ^ uJ§ {x) - ujf (x) and S{t,t-t') = 
uJq {x) J2 {x) + uJi (x), where a; = £c{t){t — t'). Thus, 
the diagonal conductivity is anisotropic with the two axis 
being defined by the linearly polarized ac-field. P{t, lo) 
and S{t,Lu) are non-zero if < 2£c{t): 



P{t,u;)^ 



TT£c{t) 

2u 
Tr£c{t) 



E(fc), 



£^c^t) 



K{k) - 2E(fc) 



(53) 



where E(fc) is the Jackobi elliptic function of the second 
kind and k = ^/T^~uJ^J4£^'(F(t) is the elliptic modulus. 
As P{t,Lo) = P{t,—uj) and S{t,iLj) — S{t,—uj) are even 
functions we work out the dc-conductivity: 



(Ttt i CT' 



vv 



P{t,Uj) 



2e'N I ^^^<! :,):'^^ S.ImG(e 

iF{e) - F(e + tj))4lmG(e + tu) - ImG(e + cj)^ 
de de 



(54) 



In the limit oj <C £c{t) the parallel to the ac-field con- 
ductivity axx is by the factor A^(t)/f ^c^(t) smaller than 
the perpendicular to ac-field conductivity ayy. This cor- 
responds to the geometry of the electron drift. In the 
absence of the electron-phonon interaction and the Lan- 
dau level mixing F{e) = v and according to Ea. H54() both 
components of the diagonal rfc-conductivity are zero. 

Electron-phonon scattering invalidates the i^(e) = v 
solution. For ep£c{t) <C i?p^, where E'^y^ = ps^ /Q'^ is the 
characteristic phonon energy, the electron-electron scat- 
tering is dominant. In this case the electron distribu- 
tion function is the Fermi-Dirac one: F{e) = l/(exp(e — 
^i)/T*) + 1), where the chemical potential /j, is deter- 
mined from the total electron normalization condition 
whereas the effective electron temperature T* is found 
from the balancing of disorder and electron-phonon scat- 
tering. A simple analysis shows that electron distribu- 
tion function is spread over many Landau levels in this 
case. On the other hand in the limit E'p^ <C e^f c(t) the 
electron-electron interaction can be neglected and weak 
electron-phonon interaction can be treated in the second 
order of the perturbation theory. In the limit V{t) ^ s 
the collision integral balancing of disorder and phonon 
scattering processes reads: 



ImG(e) j ^u(^) [F{e + w) - F(e)] ImG(e 
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IniG(e) 



s^q 



{F{e + sq) [1 - F{e)] lmG{e + sq) 



-F{e) [1 - F{e - sq)] ImG(e - sq)} 



d^q 
(2^ 



(55) 



It can be linearized by an expansion F{e) — v + SF(e): 

du) 
— 1 

277 

d^q v{l - u) s^q 



-u{uj) [SF{e + Lu)- SFie)] ImG(e + lu) 



(27r)3 



(ImG(e + sq) - ImG(e - sq)) (56) 



The r.h.s. of Ea. (|^ is evaluated as: e sQmaxi^i^ — 
iy)A{t)/Epi^, where Qmax ~ 2£c{t)/s is the cutoff mo- 
mentum corresponding to the crossover from the Lorentz 
solution to the self-consistent Born approximation solu- 
tion. Then we solve Ea. H56|l approximately as: 



5F{e) 



El^\og{{Ec{t)im) 



(57) 



This change of the electron distribution function is small 
6F{sQmax) < 1 because sQ^ax < w// < Eph. Weak- 
ness of the electron phonon interaction means that dur- 
ing the period of ac-field l/w the variation of the elec- 
tron distribution function due to the phonons is small: 
dF{t)/dt < uj. Inserting Eq.|(S7Il into Eq.(|SlI) we still 
find zero in the leading logarithmic order but in the next 
order we find approximately 

cTyy^e'Ny(l-y)uH/TEl^. (58) 

This conductivity is substantially smaller than the 'dark'- 
conductivity Oyy ^ e^N . Therefore driving electron sys- 



tem may result in that it becomes more ideal i.e. non- 
interacting with impurities system. 



VI. CONCLUSION 

In conclusion we have found that a driven disordered 
2DEG with the filling factor in high non-overlapping Lan- 
dau levels remains in a nearly coherent quantum state 
with a broad non-Fermi-Dirac electron distribution func- 
tion provided the electron-phonon interaction is suffi- 
ciently weak and the ac-electric field is strong and slow. 
In this driven state the diagonal conductivity is strongly 
anisotropic with one component decreasing with increas- 
ing amplitude of the ac-field and both being small as 
ujh /tE'^j^. The ac-field suppression of conductivity has 
a simple interpretation. Strong and slow ac-electric field 
imparts a large velocity to an electron. Fast electron ex- 
periences a reduced scattering by disorder and behave 
like a free particle. In our theory the transport scatter- 
ing time is equal to r whereas in experiments it is 100 
times larger. But qualitatively the experiment Ref^^ in 
high mobility GaAs heterostructure has found a reduc- 
tion of the dc-longitudinal conductivity in the magnetic 
fields above the cyclotron resonance loh > to that corre- 
sponds to slow ac-fields and upon increasing the power 
of microwave irradiation. 
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